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Summary. — We investigate the wealth evolution in a system of agents that ex- 
change wealth through a disordered network in presence of an additive stochastic 
Gaussian noise. We show that the resulting wealth distribution is shaped by the de- 
gree distribution of the underlying network and in particular we verify that scale free 
networks generate distributions with power-law tails in the high-income region. Nu- 
merical simulations of wealth exchanges performed on two different kind of networks 
show the inner relation between the wealth distribution and the network properties 
and confirm the agreement with a self-consistent solution. We show that empirical 
data for the income distribution in Australia are qualitatively well described by our 
theoretical predictions. 

PACS 89.65.Gh - Economics; econophysics, financial markets, business and man- 
agement . 

PACS 89.90.+n - Other topics in areas of applied and interdisciplinary physics . 



Introduction 



Empirically the literature reports several behaviors for the income and wealth dis- 
tributions in different countries. A century ago, the Italian social economist Pareto 
suggested a power-law pQ distribution in the high-income range, namely, in terms of cu- 
mulative distribution: P>(u>) oc w~ a , with a being the Pareto index On the other 
hand Montroll [3] suggested a lognormal distribution with power law tail for the USA 
personal income. More recently, wealth and income distributions in the USA and in 
the United Kingdom have been described by an exponential distribution with power law 
high-end tails @|. Whereas, the Japanese personal income distribution appears to follow 
lognormal distributions also with power law tails [3 E| ■ I n some recent papers Zipfs law 
has also been proposed y\. In this paper we add to the above empirical investigations 
an analysis of the income distribution in Australia (Figure 
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Fig. 1. - Complementary cumulative distributions for the Total annual income from all sources 
in Australia in the years 1993 — 1997. 



^From the theoretical side, it has been shown that pure multiplicative stochastic 
(MSP) processes can explain the lognormal income distribution but they fail to explain 
the power law tails [S]. Power law tails can be obtained extending MSP processes by 
including -for instance- additive noise and boundary constraints |SJ HH 1X3 ■ These 
models explain well the emergence of power law distributions, but they are incomplete, 
neglecting interactions between agents. Hence, MSPs with interacting agents connected 
through a network have been developed ^3 1141 1151 fTTTj . These models retrieve power 
law tails with exponents a which are related to the network properties. 

In this paper, we show that distributions with power law tails can emerge also from 
additive stochastic processes with interacting agents. In this case, we show that the 
network of connections among agents plays a crucial role. Indeed, the resulting wealth 
distribution is shaped directly by the degree distribution of the network. The original 
purpose of the present work was not to construct any realistic model for the wealth 
distribution. Our aim was simply to demonstrate the possibility to obtain 'fat' tails also 
without the use of multiplicative stochastic processes. Rather surprising we find out that 
the results from such an additive process are in good qualitative agreement with the 
empirical data for the income distribution in Australia. 
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Fig. 2. - Comparison between the empirical data and the theoretical (complementary cumula- 
tive) distribution associated with a scale- free network. 



2. Income distribution in Australia 

Let us briefly start with the empirical analysis of the data for the incomes in Australia. 
We analyze data from the Australian Bureau of Statistics: "Survey of Income and Hous- 
ing Costs Confidentialised Unit Record Files". In Figure Q] we report the complementary 
cumulative distributions (P>(w) = 1 — f_ p(£)d£) for the Total annual income from all 
sources in the years 1993 - 94, 1994 - 95, 1995 - 96, 1996 - 97. These data are compared 
with two possible trends in two different regions: lognormal at low and medium incomes 
and power law at high incomes. As one can see the large income region is rather well 
described with power law -like tails: P > (w) cx w~ a with exponents a respectively equal 
to 2.4, 2.6, 2.4, 2.2. Whereas the small incomes region is in better agreement with a 
lognormal distribution: P(w) = \/(ws\/Zk) exp[— \og 2 (w/x)/(2s 2 )] (with the values for 
s and x reported in the figures). Let us now introduce the theoretical framework and 
show how these behaviors can be accounted by using an additive interacting stochastic 
process. 

3. Wealth distribution from interacting additive stochastic processes 

Consider N agents which interact through a social network and suppose that at the 
time t a given agent I has a wealth wi(t). Within the same framework of other models 
proposed in the literature (141 1151 Ho] , let us first introduce a rather general expression 
for the wealth evolution: 



{l)wi(t + 1) - un(t) = A(t) + Bi(t) Wl (t) + Wj(t)Q^i{t)~ wi{t)Qi^{t) 
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where the coefficient Ai(t) is an additive noise and the factor Bi(t) is a multiplicative 
noise. These are stochastic variables which reflect market and social fluctuations. In 
addiction with these stochastic terms EquationHdescribes the exchanges between agents 
through a network: agent / receives a faction Qj_>j(t) of the wealth of agent j and 
gives a fraction Qi^j/(t) of its wealth to agent j' . The MSP model mentioned above 
takes into account only the multiplicative term Bi(t); their extensions introduce also the 
additive noise Ai (t) and the interactions Qj^i (t) . Differently, in this paper we neglect the 
multiplicative term and take into account only the additive noise and the interactions. 
In particular we assume that: i) there are no stochastic multiplicative terms = 0); 

ii) the additive term Ai(t) is a Gaussian noise with average zero and variance <7q; iii) 
each agent distributes a portion qo of its wealth equally among the other agents which 
are in contact with it through the social network. This last assumption implies: 

f 2o iff e;z> 
I elsewhere. 

where Zj is the number of agents in contact with agent j and Ij represents the set of the 
agents which exchange with agent j. Equation ^ becomes 



(3) w l (t+l) = A l (t) + (l-q )w l (t) + Y,-w J (t) . 

jeii Zj 

Note that in our case {Ai(t)) t = and Equation |3 describes a system which conserves in 
average the total wealth. 

The probability P t +\{x, l)dx that the agent I at the time t+ 1 has a wealth between x 
and x + dx is related to the probabilities to have a set {Qj^i(t)} of exchange coefficients 
and a set of additive coefficient {Ai(t)} such that a given distribution of wealth {wj(t)} 
at the time t yields, through Equation |31 to the wealth x for the agent I at time t + 1. 
This is: 

/OO /*00 f> oc 

daA t (a,l) dw\ ■ ■ j dw N P t (wj,j)S(x-a-(l-q )wi-'^2 —Wf), 
-oo J —oo J —oo -i z j' 

(4) 

where 8{x) is the Dirac delta function and At (a, I) is the probability density to have at 
time t on site I an additive coefficient Ai(t) = a. 
The Fourier transform of Equation breads: 

(5) p W fei)=^A((i-ftMi)n^Aj) ■ 

^ hi Zj 



By definition the cumulants of the wealth probability distribution are given by the 
expression: 



(6) fc M W = Hr ^_inP t( ^ Z) 



1(3=0 



where the first cumulant fc, (t) is the expectation value of the stochastic variable wi at 
the time t ((wi(t))) and the second moment k\ 2 \t) is its variance (of (£)). 
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Fig. 3. - The two degree distributions for the two networks used in the simulations. (Left: 
degree distributions, Right: complementary cumulative degree distributions). 



By taking the logarithm of Equation [5] and applying Equation [H] we get: 

(7) t\t+l)=c^+(l-g r^\t) + ^(f) V k^(t) , 

i ex, 

with cS 2 ^ = Cg and c^ v > = for any v ^ 2. This equation describes the propagation of 
the cumulants of the wealth distribution. A consequence of this equation is that the only 
moments which asymptotically can be different from zero are the first (the mean) and 
the second (the variance). With this last being directly proportional to Oq. 

We now seek for stationary solutions of Equation[7| i.e. situations in which at infinitely 
large times, the cumulants do not change in time: kj (t) = kj . 

3T. Self-Consistent solution: crystal. - Let us first consider an 'ideal' social network 
where every agent is connected with an equal number of other agents (Zj = z) in a 
perfectly ordered 'crystalline' structure. In this case, each agent is equivalent to each 
other and the asymptotic wealth distribution must be the same for every one (i.e. kj 
independent on j). ^From Equation0it follows that the expectation value for the wealth 
on each site is a constant and it is equal to the average wealth at t = 0: 

(8) («//> = fc, (1) = ^E^(°) ■ 

3 

Its variance is 

(9) of^fcW^—^S . 

l-%-(l- go ) 2 

Whereas all the other moments are equal to zero for v > 3. 

3'2. Self- Consistent solution: general case. - We now consider the more general case 
of a non-regular network. ^From Equation [7| it follows that a self-consistent solution for 
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Fig. 4. - The wealth distributions (left) and their complementary cumulative distributions 
(right) resulting from 50 simulations performed on two the different networks (1 and 2). 



the average stationary wealth on each node of the network is: 
(10) (wi) = fc ; (1) = |m , 

with m the average wealth on the ensemble of agents (m = ■ Wj (t)) and z the average 
network connectivity (z = jf z j)- Therefore the expectation value for the wealth of a 
given agent results proportional to its number of connections in the social network. On 
the other hand, we mentioned above that the only other moment which can be different 
from zero is the second. Therefore, we expect that the probability to find a given wealth 
on a given agent is a Gaussian distribution with average and finite variance. 

3'3. Wealth distribution -analytical. - The wealth probability distribution in the en- 
semble of agents is given by the sum of the distributions for each agent divided by the 
total number of agents. We have seen above that each agent has a wealth within a Gaus- 
sian distribution with average directly proportional to its connectivity (Equation I1U|) 
and finite variance. The resulting wealth distribution for the ensemble of agent is there- 
fore a weighted sum of Gaussian distributions with averages proportional to the network 
connectivity and weights given by the degree distribution. This overall distribution is 
shaped by the underlying distribution of the connectivity between agents (the degree 
distribution). It has been observed that in many social systems the degree distribution 
typically follows a power law behavior in the region of large number of connections |17j . 
This power law behavior in the social network connectivity will be therefore reflected in 
the wealth distribution which will assume a power law tail terminated by an exponen- 
tial cutoff (for a finite system). This behavior is qualitatively in agreement with what 
observed empirically. A comparison between the empirical data and the distribution 
resulting by summing a set of 3000 Gaussian distributions with averages proportional to 
the connectivity x = *xq (with xq = 30000, z — 1.28), equal variances s — 18000 and 
power law degree distribution p(z) — poz~ a (with a = 3.2, z m i n — 1, z max — 3000 ) is 
shown in Figure [21 As one can see the qualitative agreement is quite satisfactory. 

3'4. Wealth distribution - numerical simulation. - We generated large networks, with 
N = 30000 agents, by iteratively performing switching of neighbors in a triangulation 
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Fig. 5. - Comparison between the average values for the wealth of a given agent I (symbols) 
calculated from 50 simulations on the two networks (left: network 1; right: network 2) with the 
theoretical predictions from Equation 1101 (lines). 



embedded in a manifold with genus g = 10000. We introduce an 'energy' E = J2j( z j~z) 2 
and we perform a Glauber-Kawasaki type of dynamics. This procedure is an extension 
to g 7^ 1 and negative 'temperatures' of the method presented in ^S]. Two different net- 
works were generated by performing 600000 switches from a disordered start respectively 
at inverse 'temperatures' (3 = —0.5 (network 1) and j3 — +0.5 (network 2). At positive 
temperature a rather homogeneous network emerges with degree distribution centered 
around the average (z — 6 + 12(g — 1)/N) and with exponentially fast decreasing tails 
(Figure |21 network 2). On the other hand, negative temperatures favour the formation 
of inhomogeneous-scale free networks, with power law tails in the degree distribution 
(Figure network 1). 

Once the networks are generated, we associate to each agent an equal initial wealth of 
m (arbitrary) units. We set m — 100, qq — 0.1, a — 0.05/m and we run the simulation by 
updating at each time-step all the agent's wealth by using Equation|21up to a maximum 
time T. We verify that a steady state distribution is achieved after about 100 time-steps 
and therefore we set T = 1000. The resulting wealth distributions for 50 simulations over 
fixed underlying networks are reported in Figure We see that the network properties 
have a dramatic effect on the overall behavior of the wealth distribution. We observe 
an exponentially fast decay in the homogeneous network, whereas we obtain a power 
law 'fat' tail in the scale-free one. We verify that, in agreement with Equation 1101 the 
expectation velue for the wealth on each agent is proportional to its connectivity in the 
social network. Figure[5]reports the average (over the 50 simulation) wealth on each agent 
v.s. its connectivty. The theoretical prediction (Equation I10|) is also reported showing 
a remarkable agreement for both the networks. The small spreading of the data (more 
evident for network 2) indicates that non-local effects might also have some relevance. 

4. Conclusion 

We have shown that a mechanism of wealth exchange with additive Gaussian noise 
can produce distributions with power-law tails when the network which connects the 
agents is of a scale-free type. Although the original purpose of this work was not to 
produce a realistic model for the wealth evolution, we find a good qualitative agreement 
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between the empirical data and the theoretical prediction. More realistic models will 
be proposed in future works by introducing also multiplicative stochastic terms and a 
dynamical evolution in the network connectivity. 
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